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O ■ Abstract 

lO . We consider a perturbed integrable system with one frequency, and the 

approximate dynamics for the actions given by averaging over the angle. A 
classical qualitative result states that, for a perturbation of order e, the error 
of this approximation is 0(e) on a time scale 0(l/e), for e — * 0. We replace 
this with a fully quantitative estimate; in certain cases, our approach also 
gives a reliable error estimate on time scales larger than 1/e. A number 
of examples are presented; in many cases, our estimator practically coincides 
with the envelope of the rapidly oscillating distance between the actions of the 
perturbed and of the averaged systems. Fairly good results are also obtained 
in some "resonant" cases, where the angular frequency is small along the 
trajectory of the system. 

Even though our estimates are proved theoretically, their computation in spe- 
cific applications typically requires the numerical solution of a system of dif- 
ferential equations. However, the time scale for this system is smaller by a 
factor e than the time scale for the perturbed system. For this reason, com- 
putation of our estimator is faster than the direct numerical solution of the 
perturbed system; the estimator is rapidly found also in cases when the time 
scale makes impossible (within reasonable CPU times) or unreliable the direct 
solution of the perturbed system. 
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1 Introduction. 



The averaging method is a classical tool to analyse dynamical systems with fast 
angular variables: the idea is to average over the angles, to obtain an approximate 
evolution law for the slow variables (from now on, called the actions). Many appli- 
cations are physically relevant; so, error estimates for this technique on long time 
scales have an obvious interest. 

Concerning these estimates, the case of one angle is the simplest one due to the 
structure of its "resonances", which are produced only by the vanishing of the an- 
gular frequency. However, this one-frequency case covers non trivial situations: for 
example, it includes the perturbed Kepler problem, appearing in applications such 
as the dynamics of a satellite around an oblate planet and/or in presence of dragging 
(see [7] and references therein). 

The classical theory for the one- frequency case states that, under a perturbation 0(e) 
of a dynamical system with one angle and many actions, the difference between the 
actions of the perturbed and of the averaged systems is 0(e) on a time scale 0(1/ e), 
for e — > 0: see [1] [2] [3] [4] [7] (the two last references are also useful for general 
historical and bibliographical information). This is a qualitative result; the n-th 
order extensions of the averaging method proposed in the literature [4] are usually 
treated at the same qualitative level, the conclusion being that some reminder term 
is 0(e n ) on a time scale 0(l/e). To get these 0(e) or 0(e n ) bounds, one generally 
writes a number of quite rough majorizations, often containing unspecified constants 
but sufficient to obtain a linear integral inequality for the reminder; the latter is used 
to obtain the wanted bounds through the Gronwall Lemma. 

Of course, the previously mentioned results are not fully satisfactory if one aims to 
obtain precise numerical values from the error analysis; the situation is especially 
uncomfortable near resonances, i.e., when the time evolution carries the system close 
to a zero of the angular frequency. 

In this paper we show that working carefully, and avoiding unnecessary simplifica- 
tions, it is possible to derive fully quantitative and precise error estimates for the 
standard (n = 1) averaging method, for a (small) fixed e: this requires to solve 
a nonlinear integral inequality, or a related differential equation, coupled to a set 
of auxiliary differential equations. In typical cases, this is done numerically; how- 
ever, the treatment of the above system of equations is much less expensive than 
the direct numerical solution of the action- angle evolution equations; in fact, to get 
information on an interval [0, U/e) it suffices to solve the previously mentioned set 
of equations on the interval [0, U). 

To our knowledge, a quantitative error analysis for the averaging method has been 
previously proposed in [8]; however, in this reference the attention is mainly focused 
on specific applications, admitting a simple analytical treatment, rather than on a 
general scheme. In a broader sense, the present paper has some connection with [6]; 
in the cited reference, a quantitative analysis has been proposed for a rather general 



1 



class of approximation methods for the evolution equations (in abstract Banach 
spaces, so to include the case of evolutionary PDEs). 

1A. A precise setting of the problem. Let us be given an open set A of H d 
(the space of the actions) and the one- dimensional torus T: 

A = {/ = (/%!,...,<*} C R d , T := R/(2ttZ) = {tf} . (1.1) 

We fix some initial data 

Jo G A , $o £ T (1.2) 
and consider the perturbed one-frequency system 

( dl/dt = ef(l,Q) , 1(0) = J , , . 

\ dG/dt = u(l)+eg(I,e) , 0(0) = tf 1 J 

for two unknown functions I : t t— > l(t) G A, : t i— > 0(t) G T. This Cauchy 
problem contains the unperturbed frequency 

G C m (A, R) , uj{I) ^ for all / G A ; (1.4) 

the perturbation is governed by a parameter £ > 0, and by two functions 

/ = I/' ); , / e C m (A x T, R d ), 5 G C™(A x T, R), (1.5) 

(J,0).-/(J,0), <?(/,tf) ; 

throughout the paper, for technical reasons it is assumed that m ^ 2. 
From now on "the solution (I, 0) of (1.3)" means the maximal solution in the future, 
i.e., the one with the largest domain of the form [0, T), T G (0, +oo] (of course, this 
domain generally depends on the initial data). Any expression like "the solution 
(1,0) exists on D" means that D is a subset of [0, T). It is hardly the case to 
observe that I, are C m+1 functions. 

The averaged system associated to (1.3) is the Cauchy problem 

^=7(J), J(0) = /o, (1.6) 
7 = (7) G C m (A, R d ), 1^ /(/) :=-!-/ dd /(/, 0) ; 

27T J T 

the unknown is a function J : r h j(t) G A. In the same language as before, we 
stipulate that "the solution J of (1.6)" means the maximal one in the future; again, 
we have a C m+1 function. 

The system (1.6) will be compared with (1.3) for r = et, i.e., interpreting r as a 
rescaled time; if (I, 0) is the solution of (1.3) and J is the solution of (1.6) with the 
same datum I as in (1.3), the aim is to evaluate the difference 1 1— > l(t) — J(et). 
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The classical result on this subject is an estimate 

\l(t) - J(et)| < Ce for t G [0, l/e) , (1.7) 

holding for all sufficiently small e, under suitable technical conditions (especially, 
a lower bound \ou(I)\ ) c> on a convenient domain); in the above, C is a 
constant independent of e. In principle, one could obtain for C a (very complicated) 
expression, for example evaluating all the constants in the derivation of (1.7) by [1]; 
however, the explicit bound obtained in this way is not satisfactory, since in typical 
examples it largely overestimates the difference I(t) — J(et) . 

IB. Contents of the paper. Throughout the paper, the parameter e is fixed 
in (0, +oo); of course, our statements are interesting mainly if e is small (and are 
accompanied by comments which assume this). Our aim is to perform an accurate 
analysis of the distance between I and J; this will ultimately yield a bound 

\l(t) - J(et)\ ^ en(ef) for t G [0, U/e) , (1.8) 

where n : r i— > n(r) fulfils an integral inequality, or a related differential equation, 
for r within an interval [0, U). (As we will show, the existence of J, n and some more 
auxiliary functions for r G [0, U) grants the existence of the solution (I, 6) of (1.3) 
for t G [0, U/e)). 

Typically, the estimator n must be computed solving numerically the above men- 
tioned differential equation; however, this is much less expensive than the numerical 
solution of (1.3), because n depends on the "slow" time variable r = et and thus 
must be determined on an interval of length U to get an estimate for t G [0, U/e) 
(these considerations can be extended to all the auxiliary functions required in this 
approach). In the examples we will present, the function t h- > en(et) obtained 
in this way often coincides with the "envelope" of the rapidly oscillating function 
t i— > \l{t) — J(e£)|, giving practically the best possible bound of the form (1.8). Our 
bound turns out to be fairly good also in some resonant cases (where u vanishes 
at the boundary of A and the actions are close to it, either initially or over long 
times). As expected, in each example the CPU time for the computation of n is 
much shorter than the CPU time for the direct solution of (1.3). 
If U — I, Eq. (1.8) can be regarded as a quantitative formulation of the classical 
theory, involving the time scales 1 and l/e. However, in certain cases our approach 
works as well for U >> 1, yielding accurate estimates for|l(i) — J(et) \ on the extremely 
large interval [0, U/e); one can even jump to the time scales U ~ l/e, U/e ~ l/e 2 . 
The general setting of our approach is described in Section 2, where we use system- 
atically the function 

t i-> L(t) := -[I(t) - J(e*)] . (1.9) 

e 

After introducing a set of auxiliary functions and differential equations, in Lemma 
2.1 we obtain an exact integral equation for L; then, in Proposition 2.4 we derive an 
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integral inequality and show that any solution r i— > n(r) of this inequality gives a 
bound |L(t)| < n(et). For practical purposes, it is convenient to relate the integral 
inequality for n to a differential equation, which is the subject of Proposition 2.5; 
the solution n of the differential equation gives a bound \L(t)\ ^ tt(et), which is 
equivalent to Eq.(f.8). 

The subsequent Section 3 summarizes the path to n, and discusses tests for the 
efficiency of this estimator. The final Section 4 is devoted to the examples: we 
mention, in particular, the van der Pol equation, a resonant case inspired by Arnold, 
and Euler's equations for a rigid body under a damping moment linear in the angular 
velocity (which also manifest a resonance). 

To simplify our exposition, many technical aspects are treated in the Appendices. 
In particular: Appendices A, B and C contain the proofs of Lemmas 2.1, 2.3 and 
Proposition 2.5, respectively; Appendices D and E illustrate the computation of 
some auxiliary functions required by the examples of Section 4. 
The examples presented in this paper are relatively simple, since their purpose is 
mainly to test the effectiveness of the method. We postpone to later works (now 
in progress) the treatment of slightly harder applications, in particular the already 
mentioned satellite dynamics. 

2 Main results. 

2A. Some notations, i) Vectors of R d are written with upper indices: X = 
(AT l )i=i,...,(i . We use the spaces T^(R d ) of (p, g)-tensors over R d , especially for 
(p,q) = (1, 1), (2, 0) and (1,2); tensors of these three types are represented as families 
of real coefficients s/ = S3 = ¥ = (tfj k ) k = 1, d). 

Let X, Y e R d , s/, ® e T}(R d ), SB G Tjj(R d ), <€ G T*(R d ). We define the products 
XY e T 2 (R d ), s/X e R d , sf@ e T\ (R ) , <€X e T}(R d ), e R d by 

(XY) ij := X z Y j , (siX)* := <X fc , (sf®)) := , (2.1) 

(tfX)\ = tf l u X k , {^SS) 1 := ^ SS M 

(with the Einstein's summation convention over repeated indices; XX will be written 
X 2 ). We note that sf S3 is the ordinary product of si and Q) as matrices; 1^, s/^ 1 
G T}(R d ) will denote the identity matrix, and the inverse matrix of s/. The vector 
(<gX)Y = tf{XY) will be written tfXY. 

All the considered tensor spaces can be equipped with an inner product • and 
with the corresponding Euclidean norm | |. If X,Y G R d , s/,@ G T}(R d ) and 

<€,g g T^(R d ), 

d d d 

X.Y:=Y,X l Y\ s/.®=Y.^®h * (2-2) 

i=l i,j=l i,j,k=l 
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\X\ := VX.X , \si\ := sfsiTsrf , \^\ := Vtf ' •<€ . 

ii) Recalling that A C R d is open, let h : A -> R d be C". If O 1 or O 2, 
respectively, the Jacobian and the Hessian of /i at a point J are 

Let us introduce the set (open in R d x R d ). 

A t := {(/, W) G A x R d [/, 7 + 57] G A} (2.4) 

(with [7, J + 51} denoting the segment of R d with the indicated extremes. For h as 
before and £ ^ 1 or £ ^ 2, respectively, there are functions € C £_1 (A-|-, T}(R d )) 
and G C^~ 2 (A t , T*(R d )) such that 

h(I + 51) = h(I)+&(1, 51)51 , (2.5) 

ft(7 + = M') + + 1^(1, 81)61* , .**(/, 61) = Jt!«(1, 51) . (2.6) 

If d = 1, the above equations can be solved for J$? and determine them uniquely. 
If d > 1, the above equations for Jf have many solutions; in any dimension, 
explicit solutions are given by the integral formulas 

f 1 dh f 1 d 2 h 

&(I,5I):= dx—(I + x5I), J4?(I,5I) :=2 dx(l - x) — (7 + x5I) (2.7) 
Jo ol Jo "i 

(for h of polynomial or rational type, & and Jrff can be obtained more directly from 
the expression of h(I + 51)). 

In an obvious way, for a function h : A x T — > R d , we can define the derivatives 
(dh/dl)(l,#) G T{(R d ), (dh/&ff)(I,&) G R d , (d 2 h/dl 2 )(l,$) e T^(R d ) . 

iii) The average of a C l function h : A x T — > R d is the function ft, : A — > R d , 
7 I— > ft(7) := l/(27r) / t d-&h(I, •&); this notation has been already used in Eq. (1.6), 
with h — f. 

2B. The main Lemma: an integral equation for L. We consider the perturbed 
and averaged systems (1.3) (1.6), for fixed e > and initial data 7 , i?o- 
The integral equation we are going to derive will be the basic identity yielding our 
estimates on |L(t)|; it involves a number of auxiliary functions, to be introduced as 
the construction goes on. 

First of all, s G C m (A x T, R d ) and p G C m_1 (A x T, R d ) are the functions such 
that 

/ = ? + *=0; P:=|/ + | 9 . (2-8) 
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The function s, which has a preminent role in estimates on |L|, is defined by (21 
in an implicit way; an explicit formula is i 1 ) 



s = z — z 



*(/,*) :=-^ fdtf (/(/,#) -/(/)). (2.9) 

uyi) Jo 



df 

Another function to be used hereafter is the Jacobian — G C m_1 (A, T^R ')). From 
now on, U stands for an element of (0, +oo]. 

2.1 Lemma. Suppose the solution J of (1.6) exists for r G [0, U). Denote with 
R : [0, U) -> r ^ R(r) and K : [0, U) ^ R d , t ^ K(r) toe solutions of 

^ = f(j)R, R(0) = l d ; (2.10) 

^ = 5(j)K + p(j), K(0)=0 (2.11) 

(these exist and are C m ; R(t) an invertible matrix for all r G [0, [/), and K(r) = 

R(r) J^dr'R(r')- 1 p(J(0). For d = 1, R(r) = exp £ dr' %(J(r')) G (0, +oo)j. 
Furthermore, assume that the solution (I, 0) of the perturbed system (1.3) exixts for 
t G [0, C//e), wit/i (J(£*)> - J (^)) e A t- Finally, define 

L : [0 ,U/e)~B>, t»U!)~\W)-m]. ("2) 
Tnen, for t G [0, C//e) it is 

L(t) = s(l(t), 9(f)) - R(et) s(I , 7? ) - K(et)+ (2.13) 

-e(«;(I(f), ©(*)) - %(J(et)) v(I(t), ©(*))) + 

+ e 2 R(5t) £dt'R~ 1 (rf / )(«(l(t / ) ) 0(O)-|y(J(^ , ))(^ + «)(l(O ) 0(O)+ 

-^T(J(eO)u(l(0,©(0) -^(J(et'),eL(0)L(0 + ^(J(e*VL(0) L(*') 2 

In toe a&owe, v G C m (Ax T, # d ) 7 q,w £ C m ^(Ax T, R d ),ue C m ~ 2 (Ax T, R d ), and 
G C m_2 (A, I^(i2 d )) are the functions uniquely defined by the following equations: 

dv 

s = u— , v(I, -do) = for all I G A; (2.14) 



1 here, f Q means integration along any path in T from to i?; the integral depends only on 
the extremes, because the integrand has zero average. The same could be said for other integrals 
appearing later. 
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dv dv , n _ \ 

»== 5? / + w» i (2 - 15) 

dw 

p=p+u—, w(I,# o ) = for all I eh; (2.16) 
dw dw d 2 f— fdf^ 



Furthermore, $ G C m " 2 (A t , T}^)) and G C m - 2 (A t , T\(R d )) are two functions 
fulfilling Eq.s (2.5) for h = p and (2.6) for h = f : so, for (I, 51) G Af, 

p(I + 51) =p(I)+ 61)61 , (2.18) 

— — dl 1 

/(/ + W) = /(/) + -±(1)61 + -JP(1, 5I)5I 2 , se* k (1, 51) = 51) . (2.19) 

Proof. It is obtained by a long computation, where the functions s, Jif appear 
gradually. See the Appendix A. o 

Remarks, i) The above definitions must be understood in terms of the previous 
tensor notations; for example, the equivalent formulation in components of Eq. 

f) 2 ~fi fffi fffk . 

(2-i7) is ^ = dp^ijf k -^jjj- of c ° urse v(i>*) = *c,n 

w(I, tf) = cv-\I)) /* d# (p(I, #) - p(I)). 

ii) If we write l(t) = J(et) +eL(t), (2.13) becomes an integral equation for L. Most 
of the terms therein are slow, i.e., depend on et: the exceptions are L itself and 
the angle 0. The subsequent step after this Lemma will be to infer from (2.13) 
an integral inequality involving only the slow time variable et; we note that, even 
though the integral in Eq. (2.13) is multiplied by e 2 , this term appears to be of 
order e if we consider et' as the integration variable. In any case, the presence of a 
small factor e in front of the integral allows us to use for it fairly rough estimates. 

2C. A second Lemma: an integral inequality for |L|. Throughout this para- 
graph we assume that the solution J of the averaged system exists on [0, U), and 
define R, K via Eq.s (2.10), (2.11). B(I, g) denotes the open ball in H d of center I 
and radius g; we furtherly suppose the following. 

i) There is a function p G C([0, U), [0, +oo]) such that 

B(J(t), P (t)) C A for r G [0, U) . (2.20) 
We denote with T p the subgraph of p, i.e., 

T p := {(r,r)|re [0, U), r G [0, p(r)) } . (2.21) 

ii) There are functions 

a, b, c,d,eE C(T P , [0, +oo)) (2.22) 
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such that, for any r G [0, U), 5 J G B(0, p(r)) and G T, 

|s(J(r) +(JJ,t?) -R(r)s(/ ,7?o) -K(r)| ^ a(r, |<SJ|) , 



w( J(r) + 5 J, 0) - ^ ( J(r)) w ( J(r) + 5 J, 0) 



*C 6(r, |£J|) 



(2.23) 
(2.24) 



«(j(r) + 5 J, 7?) - ^y(J(r))( W + g)(J(r) + 5 J, 0)+ 



(2.25) 

(2.26) 
(2.27) 

(note that (J(r), 5 J) G Af, by the convexity of the sphere). The functions c, d, e are 
assumed to be non decreasing with respect to the second variable: 



J(t))v(J(t)+6J,&) <c(t,\SJ\) , 
\^(j(r),5J)\^d(r,\5J\) , 
\je(J(T),6J)\^e(r,\6J\) 



(r, r), (r, r') G T p , r ^ r' =>- c(r, r) < c(r, r') , (2.28) 
and similarly for d, e. Given a, 6, c, d, e, we define the functions 

a eC(T p ,[0,+oo)), a(r,r) := a{r,r) + eb(r,r) , (2.29) 

7 G C(r p x [0, +oo), [0, +oo)), 7 (r, r, £) := c(r, r) + d(r, r)£ + ^e(r, r)f . (2.30) 

We can now write the integral inequality for the function t i— > |L(£)|, with L as in 
(2.12). 

2.2 Lemma Assume that the solution (1,6) o/ the perturbed system exists on 
[0, U/e), and that \h(t)\ < p(et)/e for all t G [0, U/e). Then 

\L{t)\ ^ a{et,e\h(t)\) + e 2 \R{et)\ f dt' K^et')] 7 (et', e|L(t')|, |L(0|) . (2.31) 

Proof. We take the norm of both sides in Eq. (2.13). To estimate the right hand 
side, we use some Schwarz inequalities and Eq.s (2.23-2.27) with 5 J = l(t) — J(et) = 
eL(t); then, the thesis follows from the definitions (2.29) (2.30) of a and 7 . o 

2D. A third Lemma, on integral inequalities. To go on, we need a general 
result on a class of integral inequalities; we state it at an abstract level, forgetting 
momentarily the function |L|. 
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2.3 Lemma. Let T E (0, +00], 5 E C([0, T), [0, +00]) and 

S:={(M) I * G [0,7Ve [0,<i(t)) }, (2.32) 

H := {(t,t',£) \ t E [0,T) ,t' E [0,t], (t',£)E~}. 

Consider two functions £ E C(S, [0, +00)) and r\ E C(H, [0, +00)), £/ie Zatter non 
decreasing in the last variable: r)(t,t',£') ^ r)(t,t',£) for (t,t',£) E H and £' E [0,£]. 
Furthermore, let I E C([0, T), [0, +00)) and D G C([0, T), (0, +00)) fre snc/i t/iat 
graph I, graph D CH, and 

[(0) = , [(*) ^ £(t, l(t)) + fdt'rj^t', 1(1/)) , (2.33) 

t>(t) >f(t,t>(t))+ / dt' V (t,t',t,(t')) (2.34) 
./o 

/or aZZ* G [0, T). Taen 

I(t) < t>(f) /or aZZt G [0,T) . (2.35) 

Proof. It adapts the one of a similar result in [5]; see the Appendix B. o 

2E. The main Proposition. Throughout this paragraph we still assume that the 
solution J of the averaged system exists on [0,U), and define R, K via Eq.s (2.10) 
(2.11). We also assume there is a set of functions p, a, b, c, d, e as in paragraph 2C; 
a and 7 are defined consequently, as indicated therein. 

2.4 Proposition. Assume that there is a function n G C([0, U), (0, +00)) such 
that, for all r G [0, U), 

n(r) < p(r)/e , (2.36) 

n(r) >a(r,en(r))+e|R(r)| /" dr' {R' 1 ^ 7 (r', en(r'), n(r')) . (2.37) 

Jo 

Then, the solution (1,0) 0/ ine perturbed system exists on [0,U/e); furthermore, 
defining L as m Eg. (2.12) we have 

\L(t) I < n(et) for all t E [0, U/e) . (2.38) 

Proof. Let us recall that (1,0) is the maximal solution of (1.3), and denote its 
domain with [0, V/e); for the moment, this merely defines the coefficient V E (0, +00] 
(which can depend on e and be large, small, etc.). To go on, we provisionally put 

U' := min(V, U) ; (2.39) 



9 



one of our aims is to show that U' = U, but this will be established only in the 

second step of the proof. We also define L as in Eq. (2.12), but on the domain 

[0,U'/e). 

Step 1. One has 

\L(t)\ < n(et) for all t G [0, U'/e) . (2.40) 
To show this, we write the integral inequality (2.37) with r = et, t' = et'\ this gives 

n(et) > a(et,en(et)) +e 2 \R(et)\ [ dt 1 |R -1 (et')l 7(e*'» n ( e O) (2.41) 

Jo 

for all t G [Q,U/e), and a fortiori for t G [0, U'/e). 

On the other hand, Lemma 2.2 can be applied with the constant U therein replaced 
by U', because (1,0) is defined on [0,U'/e) and J is defined on [0, U')\ thus, Eq. 
(2.31) for |L(£)| holds for t G [0, U'/e). Now, we apply Lemma 2.3 with 

T := ^ , 5(t) := p(et)/e, (2.42) 

f (i, £) := a(et, et) , rj(t, t', £) := e 2 \R(et) \ ^(et') \ 7 (et', et, £) , 

[(*) := |L(t)| , t)(*) :=n(e*) ; 

of course, the initial condition 1(0) = holds because 1(0) = I = J(0). Lemma 2.3 
gives l(t) < D (t) , which is just the relation (2.40). 
Step 2. It is 

U' = U (2.43) 

(thus (I, G) exists on [0, U/e), and the inequality of Step 1 holds on this interval). 
It suffices to show that V ^ U; to this purpose we suppose V < U, and infer a 
contradiction. Indeed, let us put 

K := {(t,I) e [0,V/e] x R d | \I - J(et)\ ^ en(et)} . (2.44) 

This is a closed subset of R x R^; it is bounded, since t i— > J(et), t i— > n(et) are 
bounded functions on [0, V/e]. Thus, K is a compact subset of R x R d . We note 
that (*,/) G K implies / G S(J(et),en(e*)) C B(j(et), p(et)) C A (recall Eq.s (2.36) 
and (2.20)); thus, K C [0, V/e] x A. 

The previous considerations ensure compactness ofifxTcRxAxT; due to 
Step 1, we have graph (I, ©) C K x T . The inclusion into a compact set and a 
standard continuation principle for ordinary differential equations [9] imply that the 
solution (I, 0) can be extended to an interval larger than [0, V/e). This contradicts 
our maximality assumption, and concludes the proof. o 

2F. A differential reformulation of the previous results. For practical appli- 
cations, and especially for the numerical implementation of our scheme by standard 
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packages, it is convenient to replace the integral inequality (2.37) for n with a dif- 
ferential equation related to it. This equation is presented hereafter, and will be 
the basis of all applications discussed in the next sections; it is supplemented by an 
initial condition, defined implicitly by a fixed point problem. 

In the sequel we keep the assumptions at the beginning of paragraph 2E, but we 
require some more regularity on the functions a,b,c,d,e fulfilling Eq.s (2.23-2.27), 
namely, 

a,beC 2 (T p ,R) , cd.eeC^R); (2.45) 
so, the functions a, 7 in Eq.s (2.29) (2.30) are, respectively, of class C 2 and C 1 . 

2.5 Proposition, i) Assume there are real numbers £*,M ^ and a > such 
that 

£ := % - a, £* + a] C (0, p(0)/e) , (2.46) 



M < l/e , 



dot . 



^ M foreeZ , (2.47) 



\a(0,eQ -£*\ +eMa < a . (2.48) 

Then, the map £ 1— > a(0,e£) sends the interval S mto itself and is therein contractive 
with Lipschitz constant eM. So, there is a unique £ G £ solving the fixed point 
equation 

a(0,e£ )=£ . (2.49) 
£ above, let m, n G C 1 ([0, [/), i2) solve the Cauchy problem 

^ = |R- 1 | 7 (.,en,n) , m(0) = , (2.50) 

aW = V~ £ dr~ { '' £n) ) ( v 7^(-' £n ) +£ l R H R " 1 |7(-^n,n)+5|R|- 1 (R«^) ™J , 

n(0) = £ , (2.51) 

with the domain conditions 

Oct 

< n < pie , — (•, en) < l/e (2.52) 
ar 

(note that (2.50) implies m ^ 0; m £/ie above, • is the inner product of Eq. (2.2)). 

Then, the solution (1,0) of the perturbed system exists on [0, U/e) and (with L as 
in (2.12)) 

\L{t) I ^ n(et) for all t G [0, U/e) . (2.53) 

Proof. It is found in the Appendix C, after a necessary lemma. o 



11 



3 A summary of the method, and how to test it. 



3A. The main steps to implement the scheme of the previous section. In 

the approach we have outlined, the steps to be performed are the following ones. 

i) Compute J and the functions s,p, J£ , <3, M> of Eq.s (2.8) (2.14-2.19). 

ii) Determine the solution J of Eq. (1.6), on some interval [0, U); solve Eq.s (2.10) 
(2.11) for R, K on the same interval. 

iii) Find a set of functions p, a, b, c, d, e as in paragraph 2C, so as to fulfil the in- 
equalities (2.23-2.27); from them, define the functions a, 7 via Eq.s (2.29) (2.30). 
In the subsequent steps, we make on a, e the assumptions (2.45). 

iv) Determine £ , solving the fixed point problem (2.49). 

v) Search for functions m, n fulfillfing Eq.s (2.50) (2.51), with the domain conditions 
(2.52). If these equations and (1.6) have solutions on some interval [0,U), we can 
grant existence on [0, U/e) for the solution (I, 0) of (1.3), and we know that L(t) := 
(I(t) — J(et))/e fulfils on this interval the bound \L(t)\ ^ n(et). 

Here are some general comments on the practical implementation of the previous 
steps (these will also be useful to introduce the examples of the next section). 

i) Of course, the computation of /, s, p, ^ is more or less difficult depending on /, 
g and u, concerning especially the integrals over These computations can involve 
special functions (it should be noted that, in many examples coming from mechan- 
ics, /, g and u are themselves special functions). Generally, the determination of 
/, s,p, ^# is simple when, for fixed J, / and g are trigonometric polynomials in 
■d. Concerning and Jf 7 , see the remarks that follow Eq.s (2.5) (2.6). 

ii) The determination of J, R, K will be analytical in the symplest cases, and otherwise 
numerical. 

iii) For the implementation of our scheme, the functions b, c, d, e are slightly less 
important than a; in fact, they are always multiplied by the small parameter e 
whenever they appear in steps iii) iv) v). For this reason, it is important to compute 
a estimating as accurately as possible the left hand side in Eq. (2.23); as for b, e, 
in many cases one can accept rougher majorizations for the left hand sides of Eq.s 
(2.24-2.27). 

In many applications, such as in the examples of the next section, the functions 
a, b, .., e will have the form 

a(r, r) := a(j(r), R(r), K(r), r), b(r, r) : = 6(J(r), r), .., e(r, r) = e(j(r), r) (3.1) 

depending on certain known functions 

aeC 2 (A,R), 6gC 2 (T,R), c, d, e e C^f , R) , (3.2) 

with domains 

A C R d x T\(R d ) x R d x R open, T C R d x R open such that (3.3) 
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(j(r),R(r),K(r),r)GA, (J(r),r)eT for all (r,r)eT p . 
Of course, in this case it is 

a(r, r) = S( J(r), R(r), K(r), r) , 7 (r, r, £) = j{ J(r), r, £) , (3.4) 

where a e C 2 (A, R) and 7 £ C'(T x R, R) are defined by 

o?(J, St, K, r) := a(J, K, r) + eb(J, r) , (3.5) 

7 (J, r, £) := c(J, r) + d(J, r)l + ^e(J, r)£ 2 . (3.6) 

Furthermore, the derivative da/dr in Eq. (2.51) is given by 

da , . da , . dJ da , . dR da , , dK . „. 

_(, r) = _(J, R ,K,r).- + g ;f (J,R,K,r).- + g F ( J ,R,K,r).- (3.7) 

with da/dSi := (da/d&j), etc. In these situations, the function r i— > p(r) deter- 
mining the domain of a, e will often depend on r through J, i.e., p(r) = p(j(r)). 
The structure (3.1) for a, b, etc. appears naturally in cases where these functions can 
be obtained maximizing the left hand sides of Eq.s (2.23), (2.24), etc. by analytical 
means. 

In more complicated situations, one could consider the possibility to determine a, b, 
etc., maximising the left hand sides of Eq.s (2.23), (2.24), etc. by numerical (or 
partially numerical) techniques. These would give tables of numerical maxima, to 
be subsequently interpolated by elementary functions to get a, b, etc. . A second 
possibility is to derive the evolution equation for the maximum points of interest 
as function of r, to be coupled with the other differential equations in our general 
framework; this approach should work if there are no bifurcations. 
Both possibilities outlined above are especially interesting for the function a, since 
this requires the greatest accuracy; however, they will be investigated elsewhere. 

iv) The fixed point in (2.49) is given by the standard iterative formula £q = 
lim n ^ +00 / n , where l n := a(0, el n -i) and l x is chosen arbitrarily in S. One can 
compute numerically the sequence (l n ) up to a sufficiently large value n — N, and 
then assume £ ~ In- ( 2 ) 

v) Even in cases where all the other functions have known analytical expressions, 
the differential equations (2.50) (2.51) for m, n will be typically too difficult to be 
solved analytically. So, a numerical treatment will be necessary. 

If we do not have analytical expressions for J, R, K, it may be convenient to re- 
gard Eq.s (1.6) (2.10) (2.11) (2.50) (2.51) as a coupled system for the unknowns 
J, R, K, m, n, to be solved numerically on a chosen interval [0, U). 

2 By the standard theory of contractions, \Iq — In\ ^ (eM)^ 1 ^ — Zi | / ( 1 — sM), where M is 
the constant in Proposition 2.5. 
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3B. The "91-operation". Let us fix the attention on the simple situations where 
the functions /, s, J$? have known analytical expressions and a, b, c, d, e have the 
form (3.1), depending on known functions a, ...,e. It is not difficult to write a 
program of general use for these situations, which computes the fixed point £q and 
the functions J, R, K, m, n solving numerically the equations (2.49) (1.6) (2.11) 
(2.50) (2.51). From now on, the computation of £ , J, ...,n by such a program, for 
given /, ...,e (and io,#o, U), will be referred to as the Vl-operation. Of course, the 
main outcomes of this operation are the solution J of the averaged system and the 
function n binding |L(£)|. 

We have written a general program for the above purpose, using the MATHEMAT- 
ICA system. Concerning Eq.(2.51) for n, in this program the derivative da/dr is 
expressed via Eq. (3.7); the derivatives dJ/dr,dK/dr and dR/dr which occur in 
(3.7) and (2.51) are expressed via Eq.s (1.6) (2.10) (2.11) (MATHEMATICA is also 
useful, in the symbolic mode, to produce the input of the above program, i.e., the 
functions /, ...,e"; this will appear from the examples of the next section). 

3C. Testing the effectiveness of the previous method: the "£-operation". 

By the il-operation we mean, essentially, the computation of L by direct numerical 
solution of the perturbed system on [0, U/e). To avoid misunderstandings, we stress 
that in the present framework the purpose of the £-operation is merely to check the 
reliability of the estimate |L(i)| $C n(et) produced by the 91-operation, and to prove 
quantitatively that the direct solution of the perturbed system is generally much 
slower than < Jl. When U/e is very large, the ^-operation may be impossible within 
reasonable times; an example will be given in the next Section (see Figure 3f, and 
the explanations for it). Of course, the main usefulness of the 91-operation is just 
the treatment of these cases! 

To be precise, the £-operation is the numerical determination of J, L, O in the fol- 
lowing way. First, the function r G [0, U) — > J(r) is obtained solving the averaged 
system (1.6) for J; then, the functions t G [0,U/e) — > L(t),Q(t) are determined 
solving their exact evolution equations derived from (1.3) (1.6), i.e., 

f (dL/dt)(t) = f(j(et) + eL(t), Q(t)) - f(J(et)), L(0) = , ( , 

\ (de/dt)(t) = u;(J(et)) + eg(J(et) + eL(t), O(t)), 0(0) = tf . y 1 

It is easy to write a MATHEMATICA program that computes numerically J, L, O 
for given f,g,u, I ,& . 

When the ii-operation can performed within reasonable times, it can be used to test 
the 9T-procedure along these lines: 

i) one compares the graph of the estimator n (an 9t-output) with the graph of the 
function \L\ (an £- output); 

ii) one also compares the CPU times 2^, T £ for the two operations. 

These tests are presented in the next section; they are based on the programs men- 
tioned here and in paragraph 3B. In most examples, the estimator n practically 
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coincides with the envelope of the rapidly oscillating graph of |L|; furthermore, %<yi 
is generally smaller than by one or more orders of magnitude. 

4 Examples. 

In any example we consider, the initial condition for the angle is always 

tfo := . (4.1) 

Given /, g and u, the functions /, s, £f, and p, a, e are computed explicitly 
for all I (and U). After this, specific choices are made for 7 , U and e, and the 9t- 
operation is performed; to test the accuracy of the method, the ^-operation is also 
performed and some comparisons are made, as suggested at the end of the previous 
section. The results are summarized in the figures which conclude the section. Each 
figure gives the graph of the estimator n(r) provided by 91 for r e [0, U); it also 
gives the graph of |L(r/e)| in the same interval (except one case, where £ has not 
been possible within reasonable times). 

Figures referring to an example are labelled by the same number and by a letter 
(so, Fig.s la, lb and lc refer to Example 1). The legend of each figure specifies the 
choices of I , e, U, and the CPU times 1<n, (in seconds) in the execution of the 
two operations ( 3 ). 

In the chosen examples, one derives simple analytical expressions for the functions 
J, R, K but not for m, n. However, with the view of a general comparison between 
the 91- and £-operations, all examples have been treated by the general MATHE- 
MATICA programs mentioned in paragraphs 3B-3C, which solve numerically all the 
differential equations involved. Therefore, the reported times X«n, 1% include con- 
tributions from the determination of J, R, K. In any case, the analytical expressions 
of these functions are written for completeness. 
For each example: 

i) the auxiliary functions s, Jrff, p, a, e are reported in a table. All the related 
computations are analytical; the most lengthy have been performed using MATH- 
EMATICA in the symbolic mode. 

ii) The function p always gives the distance of J(r) from the boundary of the actions 
space A. 

iii) Some details on the computation of the functions a and b, c are given in the 
Appendices D and E, respectively. The expressions for d, e follow trivially from the 
ones for £f, Jj? in the corresponding tables. 

Example 1: the van der Pol equation. This is a system of the form (1.3) for 
(1,0), with 

d:=l, A:=(0,+oo), u(I) := -1 , (4.2) 
3 of course these times, depending on the PC employed, are merely indicative. 
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/(/,#) := 7(1 - ^) - Jcos(2tf) + ^cos(4tf), #(J,tf) := ^— - sin(2tf) - ^ sin(4tf) . 

The functions x := v^Icos©, v := v^IsinO fulfil the equations x = v, v = 
— x — e (x 2 — l)v, yielding the familiar van der Pol equation x + x + e (x 2 — l)x = 0. 
It is found that 

f(i) = i(i - 5) ; ( 4 -3) 

the auxiliary functions s,v, ...,J$? of paragraph 2B are reported in Table 1 ( 4 ). 
The averaged system (1.6) has the solution 

1q + (2 - 1q) e T 

for t G [0, +00), tending to 2 for r — > +00 : this long time behavior is the mani- 
festation, in the averaging approximation, of the well known limit cycle of the van 
der Pol equation (J(r) also exists for some or all r < 0, but we are not interested in 
this fact). The Cauchy problems (2.10), (2.11) for the unknown real functions R, K 
have solutions 

R(r) = (J +(9 G > K ^ = ( 4 - 5 ) 

(J + (2 - J ) e ) 

for r G [0, +00). From now on, r is confined to an interval [0, U) (and, of course, U 
will be chosen finite in the subsequent numerical computations). 
Our next step is to construct functions p, a, e as in paragraph 2C; these are also re- 
ported in Table 1 ( 5 ). All the functions (X. .... 6 clX6 C°° in (r, r), and non decreasing in 
r; they have the form a(r, r) = a(j(r), r), b{r, r) = b( J(r), r), e(r, r) = e(j(r), r), 
where a,b,c are read from the Table and d := 0, e := 1 everywhere; this corre- 
sponds to a special case of Eq. (3.1). Similar remarks could be made for the other 
Examples, but will be no longer repeated. 

Comments on this example and the figures. Figures la, lb, lc refer to the initial 
data Jo = 1/2 or Jo = 4, one below and the other above the critical value 1 = 2 (i.e., 
the limit cycle in the averaging approximation); U is 10 or 200. The ratio Tfrt/T^ 
is between 1/150 and 1/40, in the three cases. Due to the limit cycle, one expects 



4 We note that the domain Aj of ^,3^ is made of pairs (I, SI) as indicated in Table 1. In all 
the other examples, Af can be read as well from the tables. 

5 The functions 6, c constructed in this way could be replaced by appropriate, simpler majo- 
rants reducing the "confidence interval" [0, J(r)) for r; for example, one could redefine p(r) := 
min(J(r), 1/10) and infer upper bounds for b, c by means of the inequalities r k ^ r/10 k ^ 1 , for 
k = 2, 3, holding for r G [0, p(r)). These upper bounds are fairly simple, since they depend lin- 
early on r; of course, their use is correct if one checks a posteriori that < e n(r) < min(J(r), 1/10) 
for all t e [0, U). However, to perform the 01-operation in all cases presented in the figures we have 
used directly the complicated expressions in Table I, since these are easily handled by MATHE- 
MATICA. 
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Table 1. Auxiliary functions for Example 1. 



For / G (0, +oo), $ G T and 51 G (-/, +oo): 

s(I,&) = ^(4sin(2tf)-/sm(4tf)) , v{I,&) = -^8-/-8cos(2^)+/cos(4^)) , 

p(J, t?) = ^ ((4 - 2/ - J 2 ) sin(2tf) + /(/ - 4) sin(4tf) + I 2 sin(6tf)) , p(J) = , 

q(I,&) = -— (l6-10/ + 2/ 2 -(16-/ 2 )cos(2^) + /(10-2/)cos(4^)-/ 2 cos(6^)) , 
32 V / 

w(I,&) = - — (24 - 24/ - I 2 - 6(4 -21 - I 2 ) cos(2tf) + 3/(4 - /) cos(4tf) - 
96 V 

2/ 2 cos(6tf ) ) , 

u(I, #) = — — (64 - 120/ + 36/ 2 + / 3 + (-64 + 64/ + 50/ 2 - 12/ 3 ) cos(2^) + 
128 V 

+4/(14 - 17/ - I 2 ) cos(4tf) + 6/ 2 (-3 + 21) cos(6tf) + 3/ 3 cos(8tf)) , 



1/2 

J=j(r) 



128 

I 2 

= -l + /-i/ 2 , &(I,6I) = 0, Jt?(I,6I) = -1 . 

For r G [0,U), p(r) := J(r). 
For r G [0, C/) and r G [0, J(r)) : 

a(r, r ) := i ( - 2 + 10( J + r) 2 + ( J + r) 4 + 2(1 + 2( J + rff 2 ^ 

b(r, r) := — f 120 J 6 + 12/ 5 (23 + 56r) + 3J 4 (192 + 474r + 517r 2 ) + 
96 V 

+12/ 3 r(72 + 180r + 157r 2 ) + 6 J 2 r 2 (372 + 530r + 231^ + 12 Jr 3 (216 + 213r + 46r 2 ) + 
+r 4 (1404 + 690r + gir^V 72 

/J=J(r) 

c(r, r) := — (6512J 8 + 24J 7 (671 + 2096r) + 24J 6 (1693 + 5484r + 6956r 2 ) + 8/ 5 x 

x(1812 + 31188r + 39375r 2 + 38726r 3 ) + 12J 4 (768 + 4436r + 61358r 2 + 37966r 3 + 

29997r 4 ) + 8J 3 r(4680 + 39948r + 125584r 2 + 62193r 3 + 35046r 4 ) + 12J 2 r 2 (1824+ 

+52152r + 61180r 2 + 37311r 3 + 12021r 4 ) + /r 3 (119808 + 445536r + 425592r 2 + 

\ 1/2 

210995r 3 + 41976r 4 ) + 4r 4 (21600 + 33024r + 30127r 2 + 10383r 3 + 1377r 4 ) ) 

/J=J(r) 

d{r,r) : = , e(r, r) := 1 . 
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\L{r/e)\ to be bounded on the whole interval [0, +00); this fact is reproduced very 
well by our estimator n(r), that appears to approach a constant value for large r 
(see in particular Figure lc). 

Example 2: a case with action-dependent frequency. We choose 

d=l, A:=(0,+oo), u(I):=I, (4.6) 

/(/,#) := kI 2 (1 - cos(2^)) , g(I,-&) := kI 2 (1 + cos(2tf)) , k G {±1} . 
It is 

7(7) = kI 2 , (4.7) 

and the auxiliary functions s,v, ...,J%? are reported in Table 2. Let us comment on 
the vanishing of u for I — > 0. Our framework shows this "resonance" to be false: 
in fact, even though Eq.s (2.8) (2.14) (2.16) for s,v,w contain a factor 1/uj, in this 
case none of these functions is singular for I — > 0, since /, g vanish in this limit more 
rapidly than uj. 

The averaged system (1.6) is fulfilled with 

= for t e i°' w ^< w -* - { ¥h 1 1 = -1: < 4 - 8 » 

Eq.s (2.10) (2.11) for R, K have solutions 

R ( T ) - \ k(t) - < (4 9) 

on the same domain. In the sequel we assume r e [0, C/), with [/ ^ VF Kj / ; the 
functions p (the same of Example 1) and a, b, c, d, e are also reported in Table 2. 

Comments on this example and the figures. Fig.s 2a, 2b and 2c refer to the case 
k — 1, while Fig.s 2d and 2e refer to k — — 1; the initial datum is always 7 — 1- The 
two cases are radically different: in fact, according to Eq. (4.8), the solution J(r) of 
the averaged system diverges for T^l~if/« = 1, whereas for k — — 1 it is defined 
for arbitrarily large r and vanishes for r — > +00. The figures seem to indicate a 
similar behaviour for the function r 1— > |L(r/e)|; this behaviour is reproduced very 
well by our estimator n(r), which remains close to the envelope of \L(r/e) \ even for 
k — 1 and r close to 1 (see, in particular, Fig.s 2a and 2c). 
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Table 2. Auxiliary functions for Example 2. 



For I G (0, +oo), G T and 51 G (-1, +oo) : 

s(I,&) = -^/sin(2tf) , v(I,ti) = ~(1 - cos(2tf)) , 

p(J,#) = -l/ 2 (21+41 cos(2tf)+2 sin(2tf)+2J cos(4tf)-sin(4tf)) , p(I) = -i/ 3 , 
q(I, ti) = -- I 2 (2 sin(2tf) + sin(4 7?)) , 

w(J,#) = __L/^3 _4 cos (2^) +8/ sin(2tf) + cos(4tf) + 2/ sin(4#)), 
u ( J, t?) = J 2 (l6 1 2 + 10 + (40 J 2 - 15) cos(2tf) + 40/ sin(2tf) + 
+ (32 J 2 + 6) cos(4 1?) - 8 / sin(4 •&) + (8 J 2 - 1) cos(6 7?) - 8 1 sin(6 7?)' 



^(/) = 6/ 2 , &(I,6I) := -^(3I 2 + 31 51 + 5 1 2 ) , jr(J,£J):= 2k. 
For r G [0, U), p(r) := J(r). 

For t G [0,f/) and r e [0, J(r)) : 

a(r,r):=^(J(r)+r)-K(r) , 

6(r,r) := (50J 4 + (55 + 200r) J 3 + (38 + 85r + 300r 2 ) J 2 + (65 + 33r + 200r 2 ) Jr+ 
8 v 2 ^ 

\ 1/2 

+ (32 + 27r + 50r 2 )r 2 ) 

/ J=J(r) 

c(r,r) := -JL= (4608J 8 +(3904+36864r) J 7 +(1520+23296r+129024r 2 ) J 6 + (1856H- 

+5696r+ 57792r 2 + 258048r 3 ) J 5 + (4853+ 5352r+ 10032r 2 + 76160r 3 + 322560r 4 ) J 4 + 

+ (3086+ 7824r+ 11008r 2 + 56000r 3 + 258048r 4 ) J 3 r + (1862+ 2976r+ 9808r 2 + 

+21504r 3 + 129024r 4 ) J 2 r 2 + (1024+ 2312r+ 5440r 2 + 7168r 3 + 36864r 4 ) Jr 3 + 

\ 1/2 

+ (512+ 752r+ 1296r 2 + 1280r 3 + 4608r 4 )r 4 I 

/ J=J(t) 

d(r, r) := ^ (3J 2 + 3Jr + r 2 ) J=J(r) , e(r, r) := 2 . 
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Table 3. Auxiliary functions for Example 3. 



For 7 E (0, +oo), d E T and 51 E (-7, +oo) : 

s(I,#) = , v(I,ti) = — 7^(1 — cos t?) , 

I 7 2 

P( J ^) = ^2( 2sin ^ - sin ( 2 ^)) > Pi 1 ) = > 

g (J,0) = -l(3-4cos^ + cos(2^)) , w(I,#) = 
3 

u(I, ■&) = — (-10 + 15 costf - 6 cos(2tf) + cos(3tf)) , 



= , 51) := , 51) := . 

For r E [0,17), p(r) := J(r). 
For r G [0,17) and r G [0, J(r)) : 
<T,r):= ^J— - , 6(r,r):= 
d(r,r):=0, e(r,r):=0. 



(J(r)-r)3 ' 



c(r,r) := 



12 



(J(r) - r)* ' 



Example 3: a truly resonant case. Let us pass to a case where the vanishing of 
uj for 7 — > gives rise to singularities for s,v,w and other auxiliary functions. We 
assume 

d=l, A:=(0,+oo), u(I) := I , (4.10) 

/(7, i?) := 1 — cos i? , #(7,tf):=0. 

This example is considered in [4] [7] to introduce the subject of resonances; it is 
inspired by a two- frequency example in [2]. In this case, 

7(i) = i ; (4.ii) 

the functions s,....,^,J^ are reported in Table 3. The averaged system (1.6) has 
the solution 

J(r) = 7 + r (4.12) 

dl 

for r E [0, +00). Eq.s (2.10) (2.11) for R, K are very simple in this case, since — = 
and p = 0; this implies 

R(r) = 1 , K(r) = . (4.13) 

From now on, r G [0, U); the functions p, a, 6, c, d, e are reported in Table 3. 

Comments on this example and the figures. The resonance for 7 — > + could be 
expected to give problems for initial data 7 close to zero (these problems should 
appear mainly for small r, since Eq. (4.12) for J shows a departure from the res- 
onance as r grows). As a matter of fact, the estimator n approximates well the 



20 



envelope of |L(r/e)| even for small r and data fairly close to zero, such as I = 1/2 : 
the agreement is rather good for e = 10~ 2 (Fig.s 3a and 3b) and very good for 
e = 1(T 3 (Fig.s 3c and 3d). 

The agreement between n and the envelope of |L| is very good even for e = 10 -2 , if 
we consider the larger datum I — 2 (Fig.3e). Fig.3f refers to the same situation on 
the larger interval r G [0,200). The statement on X £ in the legend means that the 
numerical computation of L was interrupted after 240 seconds, when the package 
had not yet produced a result; note that, on the contrary, the ^-operation for the 
same interval is very fast. 

Example 4: damped Euler's top. We consider the system (1.3), with 

d = 2, A := {I = (J 1 , J 2 ) | I 1 , I 2 G (0, +oo)}, = I 1 1 2 , (4.14) 

/(/,#) := (-/ 1 (A 1 + /iCOs(2^)),-/ 2 (A 2 -/icos(2^)) , g(I,#) :=/isin(2tf) ; 
this depends on three real coefficients /i, Ai, A 2 for which we assume 

Ai > , -Ai < /i < Ai , A 2 > -Ai . (4.15) 

This system is related to Euler's equations for the components p, q, r of the angular 
velocity of an axially symmetric top, in presence of weak damping. More precisely, 
assume that the moment of the damping forces is a linear function of the angular 
velocity, and that the linear operator expressing this dependence has a diagonal 
matrix — e diag(£', F, G) in the reference system in which the inertia operator has 
the form diag(A, A, C), with A, C, E, F,G,e G (0, +oo) ( 6 ). Then, Euler's equations 
are 

Ap + (C-A)qt= -eEp , Aq - (C - A)pt = -eFq , Cx = -eGx . (4.16) 

Given this system, we define //, Ai, A 2 through the equations 

E = A(fi + \ 1 ), F = A(\ 1 - I i), G = C(Ai + A 2 ) , (4.17) 

which imply the inequalities (4.15). Now, if (1,0) = (I 1 ,! 2 ,©) is such that I = 
0) and 6 = uj(I) + sg(l, 6), the functions 

A 

p^I^cosO, q^I^inG, t := -I 1 ! 2 (4.18) 

O — A 

fulfil Euler's equations (4.16). 

Let us return to (4.14). This implies 

TOO = ("Ai/\ -A 2 / 2 ) ; (4.19) 

6 Of course, quantities like A, G, the time t, etc., can be treated as real numbers, because we 
suppose to have fixed all the necessary physical units. 
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the functions s, Jt? are reported in Table 4. The averaged system has the 
solution 

J^r) = P e~ x ' T (i = l,2) (4.20) 

for r e [0, +oo). Eq.s (2.10) (2.11) for the 2x2 matrix function R and for the 
2-component function K have the solutions 

R(t) = diag(e- AlT , e" A2T ) , K(r) = (0,0) . (4.21) 

From now on, r is confined as usually to an interval [0, U). The functions p, a, e 
for this example are reported in Table 4; the length of the expressions of b, c is 
mainly due to the need for covering all possible values of Ai, A2,//. 

Comments on this example and the figures. In this case the main difficulty is the 
fact, following from (4.20) (4.15), that J : (r) J 2 (r) = l£ I 2 e~^ 1+ ^ T is small for large 
t. On the other hand, oj(I) vanishes for I 1 I 2 — > 0, and in this limit many auxiliary 
functions diverge; so, the averaged system falls exponentially into a resonance. 
In this situation one expects a rapid growth of |L|, which is in fact confirmed by 
Fig.s 4a-4d; the same figures show that our estimator n(r) approximates well the 
envelope of \L(t/e)\ on [0, U], when U is of the order of the unity. In Fig. 4d, a good 
agreement between |L(r/e)| and |n(r)| is attained on the longest interval among the 
four pictures (namely, for r G [0,3)). This is because we take, simultaneously, the 
largest value for 7q Iq and the lowest values for e and Ai + A2. 



Table 4. Auxiliary functions for Example 4. 



For / = (J 1 , 1 2 ) E (0, +00) 2 , ■& E T and 51 = (SI 1 , 5I 2 ) E (-J 1 , +00) x (-J 2 , +00) : 
s(I, = £ sin(2^) (-1 , 1) , v(I, *) = sin 2 (-1 1 

P (I, = ^1 j- ^M , Al + 3 7 S ^ ) , p(J) = (0, 0) 
_/xsin 2 tf ^ 2A 2 + 2/i + Ai + /i cos(2^) A 2 + 2p + 2A X + 3p cos(2^) ' 

9v , 17 J 



2/1/2 I p ' /I 
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IT n . usin 2 ^ / \ 2 + ucos 2r d X 1 + 3a cos 2 #\ . T n . usin 2 ^ 

W ( J ^) = -^fTp- [ p > 71 J > m ( / ^) = 7j t 72- x 

/ 4A 2 . + 6A 2 /i + 2A 2 Ai + /iAi + /i(4A 2 + 3/x + Ai) cos(2$) + 3^ 2 cos 2 (2$) 

X l P ' 

3A 2/ u + 2A 2 Ai + 10/iAi + 4A 2 + 3^(A 2 + 5/x + 4A X ) cos(2tf) + 15/i 2 cos 2 (2tf) \ 

T ^ / ' 

f (/) = ( " Al ^ ) , Jt{l) = ( " A ? _° A2 ) , Sf(J,tfi) = 0, ^(1,81) = 

For t G [0, 17), p(r) := min(j 1 (r), J 2 (r)). 
For r G [0, [/) and r G [0, J(r)) : 

1/2 



o(r,r) := — ( 7TT73 — + 

b(r,r) : 



2 \(3 l (r)-r) 2 (J 2 (r)-r) 2 

(ftnJ^r) 2 + 6 22 J 2 (r) 2 + b^\r)r + 6 2 J 2 (r)r + b P) 1/2 



8(J 1 (r) -r) 2 (j 2 (r) - r) 2 
fen := 16(A 2 + A 2 ) + Ai(12A 2 + 20|A 2 |) + 2(Ai + A 2 )/x + 4(Ai + |A 2 |)M + /i 2 
6 22 := 16(A 2 + A 2 ) + Ai(12A 2 + 20|A 2 |) + 6(A X + A 2 )/i + 12(A X + |A 2 |)|/x| + 9/x 
61 
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32(A 2 + A 2 ) + 64A!|A 2 | + 12{X 1 + |A 2 |)|/x| + 2/i 2 , 
32(A 2 + A 2 ) + 64A!|A 2 | + 36(Ai + |A 2 |)|/i| + 18/i 2 , 
16(A 2 + A 2 ) + Ai(12A 2 + 20|A 2 |) + 4(Ai + A 2 )/i + 14(Ai + |A 2 |)|//| + 9// 2 

(cnJ^r) 2 + c 22 J 2 (r) 2 + dJ^Tjr + c 2 J 2 (r)r + c r 2 ) 1/2 



C(T ' r):_l/i| 32(j 1 (r)-r) 2 (j 2 (r)-r) 2 
c n := 1024(A| + A 4 ) + 6144A 2 A 2 + 512(A 2 + A 2 )A!(3A 2 + 5|A 2 |) + 640(A 3 + \ 3 2 )n+ 
+896(A?+|A 2 | 3 )|/i| + 1920(Ai + A 2 )A 1 A 2 /i + 2688(Ai + |A 2 |)Ai|A 2 ||/i| + 704(A 2 + A 2 )/i 2 + 
+32A 1 (17A 2 + 27|A 2 |)// 2 - 24(Ai + A 2 )/i 3 + 264(Ai + |A 2 |)H 3 + 27/i 4 , 
c 22 := 1024(At + A 2 ) + 6144A 2 A 2 + 512(A 2 + A 2 )A!(3A 2 + 5|A 2 |) + 384(A 3 + A 3 )/i+ 
+1408(A 3 +|A 2 | 3 )|/i| + 1152(A 1 +A 2 )A 1 A 2 /i+4224(A 1 +|A 2 |)A 1 |A 2 ||/i|+2816(A 2 +A 2 )/i 2 + 
+32A 1 (21A 2 + 155|A 2 |)// 2 + 120(Ai + A 2 )/i 3 + 1800(Ai + |A 2 |)H 3 + 675// , 
ci := 2048(A 4 + A 4 ) + 12288A 2 A 2 + 8192(A 2 + Al)A 1 |A 2 | + 3072(A 3 + |A 2 | 3 )|/i| + 
+9216(A 1 + |A 2 |)A 1 |A 2 || / u| + 1408(A 2 +A 2 ) / u 2 +2816A 1 |A 2 | / u 2 + 576(A 1 + |A 2 |)|/i| 3 +54/i 4 , 
c 2 := 2048(A 4 + A 4 ) + 12288A 2 A 2 + 8192(A 2 + A 2 )A!|A 2 | + 3584(A 3 + |A 2 | 3 )|/i| + 
+10752(Ai + |A 2 |)Ai|A 2 ||/i|+5632(A 2 + A 2 )/i 2 + 11264Ai|A 2 | / u 2 + 3840(Ai + |A 2 |)|/i| 3 + 
1350// , 

c := 1024(A 4 + A 4 ) + 6144A 2 A 2 + 512(A 2 + A 2 )A!(3A 2 + 5|A 2 |) + 512(A 3 + A 3 )/i+ 
+2048(A 3 +|A 2 | 3 )|/i| + 1536(A 1 +A 2 )AiA 2 /i+6144(Ai + |A 2 |)A 1 |A 2 ||/i|+2816(A 2 +A 2 )/i 2 + 
+32Ai(19A 2 + 157|A 2 |)/i 2 + 48(Ai + A 2 )/i 3 + 1872(Ai + |A 2 |)|/i| 3 + 675// ; 

d(r,r) :=0 ; e(r,r):=0. 
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Figure la. I = 1/2, e = 1(T 2 , U = 10. 
Graphs of n(r) and |L(r/e)|. = 0.062s, 
T £ = 3.2s. 
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Figure lc. I = 4, e = 10" 2 , C/ 
Graphs of n(r) and |L(r/e)|. X«n 
% z = 67s. 
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0.1 0.2 0.3 0.4 0.5 

Figure 2b. k = l,/ = 1,£ = 10~ 2 (as in 
Fig. 2a). Graphs of n(r), |L(r/e)| in a de- 
tailed view, for r G [0, 0.5]. 




Figure lb. I = A, e = 10 

Graphs of n(r) and |L(r/e)|. 
T £ = 3.0s. 



10 
= 10. 

0.078s, 




0.6 

Figure 2a. « = 1,J = M = 10~ 2 ,U = 0.9. 
Graphs of n(r), |L(r/e)| (note that J(t) — > 
+ooforr^l-). 1^ = 0.0325, T £ = 0.36s. 
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0.76 0.78 0.82 0.84 0.86 0.88 0.9 

Figure 2c. k = 1,Iq = l,e = 10 -2 (as in 
Fig. 2a). Graphs of n(r), \L(r/e)\ in a de- 
tailed view, for r G [0.75,0.9]. 
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Figure 2d. k = -1,I = M= 10~ 2 ,f/ = 200. 
Icn = 0.078s, T £ = 0.58s. Graphs of n(r), 
|L(r/e)| in a detailed view, for r e [0, 10]. 
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Figure 3a. I = 1/2, e = 10~ 2 , C/ = 10. 
Graphs of n(r) and \L(r/e)\. 1<yx = 0.23s, 
% z = 0.95s. 
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Figure 3c. /„ = 1/2, e = 10" 3 , C/ = 10. 
Graphs of n(r) and |L(r/e)|. X<n = 0.23s, 
X £ = 12s. 




Figure 2e. /c = — 1,J = l,e = 10 2 (as 
in Fig. 2d). Graphs of n(r), |L(r/e)| in a 
detailed view, for r G [100,200]. 
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Figure 3b. I = 1/2, e = 10~ 2 (as in 
Fig. 3a). Graphs of n(r) and |L(r/e)| in a 
detailed view, for r e [0, 1). 
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Figure 3d. I = 1/2, e = 10~ 3 (as in 
Fig. 3c). Graphs of n(r) and |L(r/e)| in a 
detailed view, for r e [0, 1). 
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Figure 3e. I Q = 2, e = 10~ 2 , U = 10. 
Graphs of n(r) and |L(r/e)|. 1<y\ = 0.16s, 
T £ = 1.2s. 
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Figure 4a. n = 1, A x = 2, A 2 = -1, = 4, 

J 2 = 4, e = 10~ 2 , Z7=l. Graphs of n(r) and 
|L(r/e)|. 0.047s, X £ = 1.7s. 
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Figure 4c. fi — 1, Ai = 2, A 2 = -1, — 4 , 

/2 = 1^ = 10-3, [/■=!. Graphs of n(r) and 
|L(r/e)|. X<n = 0.047s, T £ = 4.2s. 
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Figure 3f. I = 2, e = 10" 2 , U = 200. 
Graph of n(r). 1^ = 0.28s, X £ > 240s. 
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Figure 4b. p = 1, Ai = 2, A 2 = -1, Iq 1 = 4, 
7 2 = 1,£ = 10- 2 ,C/=1. Graphs of n(r) and 
\L(r/e)\. Z m = 0.047s, X £ = 0.44s. 
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Figure 4d. // = 1, Ai = 1.1, A 2 = -1, 1% = 4, 

J 2 = 4, e = 10~ 3 , C/ = 3. Graphs of n(r) and 
|L(r/e)|. Z m = 0.046s, T £ = 99s. 
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A Appendix. Proof of Lemma 2.1. 

First of all, the Cauchy problem (2.10) has a (unique) solution on [0, U), and this 
is C m , because we have a linear differential equation for R, with a C m ~ x matrix 
df 

function r i— > — (J(t)). The invertibility of R(r) follows from the Wronskian identity 

df 

detR(r) = detR(O) exp dr' tr — (J(r')) and from the initial condition R(0) = 1^; 

the d — 1 expression of R is obvious. The statements on K that follow Eq. (2.11) are 
also elementary (as for the C m regularity, note that p(J) is a C m ~ x function of r). 
To go on, we introduce the short-hand notations 

J,R, K, ^,etc. = the functions 1 1— > J(et) , R(et) , K(et) , ^-(et) , etc. ; (A.l) 

(XT (XT 

in the same spirit, for /i:AxT^ R d and fc : A — > R d we also intend 

/i, fc, fc(J) = the functions i i-> 9(t)) , 1 1-> fc(l(t)) , t i-> fc(J(et)) . (A.2) 

In these notations, one has L = (I — j)/e and Eq.s (1.3) (1.6) imply 

we continue dividing the argument in steps. 
Step 1. One has 

dL ds df , . 1 9 ... , o 

*= w w + e a7 (J)L+ 2 e ^ (J ' EL)L • (A4) 

In fact, Eq. (A.3) and the first Eq. (2.8) imply 

dL (9s — — 

now, it suffices to recall that I = J + eL and use Eq. (2.19) with (I, 51) replaced by 
(J,eL) . 

Step 2. For each function heC\Ax T, R d ), it is 

dh dh fdh , dh \ , . ^ x 

u w = it- e \m f+ m 9 )- (A - 6) 

This follows easily from 

dh dh dl dh dQ dh , dh , . , . 

Tt = dfdt + Wit = £ di f + M {UJ + £g) ■ (A ' 7) 
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Step 3. One has 

ds ds dw 

To prove this, we note that Eq. (A. 6) with h — s and the definition (2.15) of p give 

ds ds , . v 

^ = Tt- cp - (A ' 9) 

On the other hand, Eq.s (2.16) (A. 6) with h = w and the definition (2.17) of u imply 

_ dw _ dw , . i n\ 

p = p + UJ ~d~d =p + Itt ~ £U] (A.10) 

furthermore, Eq. (2.18) with (1,51) replaced by (J, eL) gives 

p = p(j)+e^(J,eL)L. (A.ll) 

Inserting Eq. (A.ll) into (A. 10), and the result into (A. 9), we get the equality (A. 8). 
Step 4. One has 

ft ~ £ § {J)L = ft- £ ^dJ- £p{J) + £2{u ~ ^ £L)L + \^ { J) £L)l2) • (A ' 12) 

This follows immediately from Eq. (A. 4) and from the equality (A. 8). 

Introducing the next steps. Eq. (A. 12) is an equality involving total derivatives, 
and nonderivative terms proportional to e or e 2 . Our aim is to obtain an equality 
for L involving only total derivatives and nonderivative terms proportional to e 2 ; 
due to the structure of the terms in e of Eq. (A. 12), this result can be achieved 
using the functions R and K. In the sequel we will derive some identities involving 
R, where the operator R(rf/rft)R _1 plays a major role; inserting these relations into 
Eq. (A. 12) (and factoring out R) we will finally obtain an identity with the desired 
structure, where the nonderivative terms are confined to the order e 2 . 

Step 5. One has 

^ - « ■ < A13 » 

For each C 1 function X : [0, U/e) -> R d ,t^ X(t), this implies 

Eq. (A. 13) follows from the relation 

d. u dR , dR' 1 df , x dR- 1 . A 

= * (RR > = * R + R ^T = E a7 (J > + k ~m • < A15) 
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where, in the last passage, we have used Eq. (2.10) to express dR/dt = sdR/dr. 
Having extablished (A. 13), we consider any function X as above and note that 

dX df , x dX dR- 1 , ,dX dR' 1 . 

— - e-V J X = — + R^— X = R(R — + — — X , A.16 

dt dr 1 dt dt v dt dt ' v 1 

whence Eq. (A. 14). 
Step 6. One has 

dh d f . s d . i s dw d , i . idf,, , . „ „, 

* - e a7 (J)L = V R L) • e * = eR S (R w) + E 57 (J > ■ (A17) 

ep(J) = R^-(R _1 K) (A.18) 



dt 



(A.19) 



(note thai tie rigM nand sides of Eq.s (A.17-A.19) all appear in Eq. (A. 12)). 
Eq.s (A. 17) are mere applications of the general identity (A. 14) with X = L and 
X = w (i.e., the function u>(l,©)), respectively. Eq. (A.18) follows writing (A. 14) 
with X = K, and expressing dK/dt = edK/dr via Eq. (2.11) . The derivation of Eq. 
(A.19) is a bit longer. First of all, from Eq. (A. 14) with X = s we infer 

ds d . i . df . , / . s 

*= R S (R s) + E 57 (J)si < A20) 

to continue, we will reexpress — — (J)s as Rx a total derivative, plus terms of the first 

dl 

order in e. To this purpose, we write s in terms of v via Eq. (2.14), and then use 
Eq. (A. 6) with h — v; this gives 

dv dv /dv „ dv \ dv ,. _ x 

• = u « = *- e (fl7 /+ « 9 ) = *- e9 ' < A21) 

the last passage following from the definition (2.15) of q. This implies 

f <^ - - 1 (f < j >) - ■ < a - 22 > 

dj 

On the other hand, Eq. (A. 14) with X = ~qj(^) v an d Eq. (1.6) give, respectively, 
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s(w (J) ) = 9F (J) A =E aF (J)/(J) - (A - 24) 

Substituting Eq.s (A.23-A.24) into (A.22), and recalling the definition (2.17) of Jt , 
we finally get 

M(J) S = £ 4 f R_1 |?( J ) A - e 2 (j*{Z) v + 2( J) q] ; (A.25) 



or ' dt\ dr ' ) V «9/ 

inserting this result into Eq. (A. 20), we obtain the desired relation (A. 19). 
Step 7. One has 

± (R-L) = ± (E-'(. - K)) - 4< R- >(„ - |( JH ) + (A.26) 
97 l 

+e 2 R _1 ( u - ^y(J)(™ + g) - ^#(J)v - Sf(J, eL)L + -J^(J, eL)L 2 ) . 

To prove this, we return to (A. 12) and reexpress dL/dt — e— (J)L, ds/dt, edw/dt, 

ep(J) via Eq.s (A.17-A.18). Multiplying both sides by R" 1 , we obtain Eq. (A.26). 
Step 8. Conclusion of the proof. We integrate Eq. (A.26) from to t, explicitating 
the dependence of all objects on I, 0, J, t and taking into account the initial condi- 
tions for I, O, J, R, K, as well as the relations L(0) = 0, v(I, i?o) — i?o) — 0. This 
gives an expression for R _1 (et) L(t): multiplying by R(et), we get the thesis (2.13). 



B Appendix. Proof of Lemma 2.3. 

As anticipated, we are inspired by the proof of a similar statement in [5] (see Chapter 
XII, § 23, Theorem 1); therefore we merely sketch the argument. Let us define 

& ■ = {t x e (0,T) | [(*) < t)(t) for all t E [0,*i) } , Ti := sup & . (B.l) 

(Note that (2.34) and (2.33) give t>(0) > f (0, t)(0)) > = [(0); so, by continuity, ,9 
is nonempty). In the sequel we will assume T\ < T, and infer a contradiction. 
From (B.l), it is clear that [(Ti) < t)(Ti). We cannot have [(Ti) < t)(Ti) since this, 
by continuity, would be against the definition of Ti; thus 

[(Ti) = t)(Ti) . (B.2) 

On the other hand, the assumptions of the Lemma and (B.2) imply 

k t i) m, kti))+ [ 1 dt'v^t', [(0) = ( 2) em, t»(Ti))+ / 1 ^(T^f, [(0) ^ 

Jo Jo 
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< (3) £(TiMTi)) + / dt' V (T u t',»(t')) < (4) 0(71) , (B.3) 
Jo 

which gives again a contradiction. (For better clarity: the relations (1)(2)(3)(4) 
follow, respectively, from (2.33), (B.2), the monotonicity of rj and (2.34)). 



C Appendix. Proof of Proposition 2.5. 

We begin with a Lemma; this holds under the same assumptions written at the 
beginning of paragraph 2E, before stating Proposition 2.4. 

C.l Lemma. Assume that there is a family of functions n$ G C([0, Us), (0, +oo)) ; 
labelled by a parameter 5 G (0, 5J, such that the following holds: 

i) U & ^U for 5 -> 0+; 

ii) for all 5 G (0, 5J and r G [0, Us), it is 

n s (r) < p(r)/e , (C.l) 

n s (r) =5 + a(T,ens(T))+s\R(r)\ f dr' K^r')] 7 (r', exis(r'), n s (r')) ; (C.2) 

Jo 

Hi) for each fixed r G [0, U), the limit n(r) := lim < 5_ +0 + n s( r ) exists in [0, +oo) (note 

that t G [0, Us) for sufficiently small S, due to i)). 

Then the solution (1,0) of the perturbed system exists on [0,U/e) and 

\L(t)\ <: n(et) for all t G [0, U/e). (C.3) 

Proof. Of course, ii) implies 

n s (r) >a(r,Ens(r))+e\R(r)\ f dr' \*t\t')\ 7 (r', en s (r f ), xi s (t')) (C.4) 

Jo 

for all 5 G (0,5*] and r G [0,Us). Therefore, Proposition 2.4 can be applied to the 
function n<5 on the interval [0, Us); this implies that (I, O) exists on [0, Us/e), and 

|L(t)| < n s (et) for all t G [0, U s /e). (C.5) 

Now, sending 5 to zero and using hi) we easily obtain the thesis. o 

We now pass to Proposition 2.5. So, we have the assumptions at the beginning of 
paragraph 2E, strengthened by the smoothness requirements (2.45) for a, b, c,d,e . 

Proof of Proposition 2.5. For the sake of brevity, we put 

a : £ -> R, t ^ ao(t) ■= «(0, el) (C.6) 
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and extend this definition to any 5^0 setting 

a s : E -> R , ^ a 5 W := a (^) + 5 • (C.7) 
We proceed in several steps. 

Step 1. For each 5 ^ 0, as is a contractive map. In fact, for all £, £' G S, we have 

|aj(*)-a*(OI=e 

But eM < 1 by the first inequality (2.47), so the thesis is proved. 

Step 2. There is 5* > such that, for all 5 G [0, 5*], sends E into itseif. In fact, 

for any 5^0 and £ G E, 

M*) - t| = KW + 5 - L\ < \a (£) - a (Q\ + \a (Q - £*\ + 6 ^ 

^ eM\£ - Q + \a (Q - 4| + S ^ eMa + \a (Q - Q + 5 , (C.9) 
where the second inequality follows from Eq. (C.8) with 5 = 0. Now, let us define 

5* := (1 -eM)a- \a (Q - Q , (CIO) 

and note that 5, > by (2.48). For 5 G [0,5*] and £ G E, Eq.s (C.9), (CIO) imply 

- £*\ < o", i-e., a*(^) e e. 
Step 3. For all 5 G [0,5*], tie map ct^ has a unique fixed point £5 G E, which 
depends continuously on 5. Existence and uniqueness of the fixed point follows 
from the Banach theorem on contractions; to prove continuity we note that, for all 
5, 5' G [0, 5*], 

[4 - 4'| = \a s {U) - = |a (4) + ^ - «o(40 - 5'| «S (C.ll) 

^ \a (£ 5 ) - a (£ 5 >)\ + [5 - 5'\ < eM\£ s - £ 5 >\ + \6 - 5'| , 
the last inequality depending on (C.8) with 5 = 0. This implies 

(C.12) 

so the map 5 ^ £s is Lipschitz, and a fortiori continuous. 
Step 4. Proving the thesis of i). This follows from Step 3, with 5 = 0. 
Step 5. Proving the thesis of ii). For any 5 G [0, 5*], let £s be as in Step 3. From the 
standard continuity theorems for the solutions of a parameter-dependent Cauchy 
problem, we know that there is a family (Us, tris, ti$)<5e(o,<s,] with the following prop- 
erties a) b): 



da . 



\£-£'\ < eM|£-£'| 



(C.8) 



32 



a) for all 5 G (0,5*}, it is m^n* G C 1 ([0, Us), R); furthermore, these functions fulfil 
the equations 

= |R- 1 | 7 (-,£n 5 ,n (5 ) , m 5 (0)=0, (C.13) 



dr 

dris /„ <9a . .x- 1 f da . . . .. _■, . . . . , dR, \ 

~dr = \ ~dr ('' en %) ( (•> en «) + e l R ll R I 70, £n *> n*) + e|R| (R.— jm* I , 

na(0) = 4 (C.14) 

with the domain conditions 

do. 

0<n 5 <p/e, —(-,ens)<l/e. (C.15) 
or 

b) One has 

U s ^ 5 ^o+ U, n s (t) ^^ + n(t), m s (t) -^ + m(f) for all i G [0, 17), (C.16) 
where m, n are as stated in ii). 

Let us consider the pair m.,5, for any 5 G (0, 5*]. Then, integrating (C.13), 

m*(r)= rdr , |R- 1 (T , )|7(r',£n 5 (r , ),n 5 (r / )) for r G [0, 17*). (C.17) 
Jo 

Furthermore, from Eq. (C.14) we infer 

I ' (9a. N \c?ri5 /(9a. . . .. _ 1( . , . , dR, \ 

0= \l-e—(-,ens))—-[-^(-,ens)+e\R\\R \ IV, en*, n*) + e|R| ( R, ^l) m <5 J 

/ 9a. ,\dm /da. s . .drru d\R\ \ /^,„„s 

= { 1 -^^ ) H-W^ + ^ +£ ^) • (ai8) 

the last passage depends on Eq. (C.13) for m«5, and from the identity d\R\/dr = 
dy/R»R/dr = |R| _1 (R.rfR/rfr). The result (C.18) can be rephrased as 

= ^-(n 5 - a(-,en 5 ) - e\R\m s ) ; (C.19) 
dr 

the constant value of the above function can be computed setting r = 0, and is 

n 5 (0) - a(0, en*(0)) =i s - a(0, ei s ) = U - «o(4) = S (C.20) 

(recall the initial condition in (C.14), Eq.s (C.6) (C.7) and Step 3, giving £ s = 
a s{^s) — Oio(£s) +S). Therefore, 

ns(r)-a(T,ens(r))-e\R(T)\ms(r) = 5 for r G [0, Us) ■ (C.21) 

From Eq.s (C.21) and (C.17), we see that 1x5 fulfils Eq. (C.2) of Lemma C.l. Due to 
Eq. (C.16) on the limit for 5 — > + , from Lemma C.l we finally obtain the thesis, o 
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D Appendix. The functions a of the examples. 



Example 1. One must determine a function fulfilling Eq. (2.23) for r G [0,U), 
5 J G (— J(t),J(t)) and i? G T. Neither K nor R (nor the initial datum) play a 
significant role in this computation, since K = 0, s(/o,$o) = and R(r) appears in 
Eq. (2.23) as a multiplier for the second of these vanishing terms. In conclusion, to 
obtain a we can simply bind |s(J(r) + SJ,"d)\ in terms of J(r) and r := \SJ\. 
Consider any point 7 G A; of course, 

max|s(7,tf)| = (maxs 2 (7,tf)J . (D.l) 

Derivating with respect to i?, one finds that the maximum of s 2 is attained for 
cos 2 $ = 1/2 + (1 — \/l + 27 2 )/(47); by elementary computations, this gives 

max|s(7,tf)| = a(7) , (D.2) 



where a is the C°°, strictly increasing function given by 

a : (0, +oo) -> (0, +oo) , / h-> a(J) := ^ ( - 2 + 10/ 2 + J 4 + 2(1 + 27 2 ) 3 / 2 ) 1/2 . (D.3) 

Let r G [0, [/), G (-J(r), J(r)), 7? G T and r := |<$J|. Then, 

|s(J(r) + (JJ,7?)| ^a(j(r) + 5J) ^ a( J(r) + r) ; (D.4) 
the last term above is just the function a(r,r) of Table 1. 

Example 2. We refer again to Eq. (2.23); as in the previous example, R(r) plays 
no role, since it appears in Eq. (2.24) as a multiplier for the term s(/ , #0) = 0. A 
simple computation gives 

\s(J(t) + 6J,&) -K(r)| = I -£(J(r) + <JJ)sm(2#) -K(t)| ^ (D.5) 

< ^(J(r) + + |K(r)| = i(J(r) + |5J|) " K(r) ; 

this means that Eq. (2.23) is fulfilled by the function a in the Table 2. 

Example 3. Again, R and K play no role in the analysis of Eq. (2.23) and it suffices 

to bind |s(J(t) + 8 J, i?)| in terms of r := \SJ\. Clearly, 

therefore, Eq. (2.23) is fulfilled by the function of Table 3. 

Example 4. In the left hand side of Eq. (2.23), the terms K(r) and s^o,^) are 
zero; so, to find a we must bind |s(j(r) +5J, in terms of r := \5J\. Maximization 
with respect to d can be done analytically; as a final result, Eq. (2.23) is fulfilled 
by the function a in Table 4. 

In each example, the function a determined as above gives an accurate bound on 
the left hand side of Eq. (2.23). 
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E Appendix. The functions 6, c of the examples. 



i) In comparison with a, the functions b, c, d, e in Eq.s (2.24-2.27) can be constructed 
using rougher majorizations, see the comments in paragraph 3A. Here we fix the 
attention on b and c, since the functions d, e of the examples are obtained trivially. 

ii) In all the examples, to find b and c we must essentially derive a majorant for an 
expressions of the form h(J(r), 5 J, where h(J, 5 J, i?) is a trigonometric polynomial 
in whose coefficients are polynomials in J and 5 J. The majorant should depend 
only on J(r) and \8J\; so, the problem is reduced to finding a function k such that 

h(J, 5 J, #) ^ k(J, r) for e T and r := |5J| . (E.l) 

Let us exemplify this situation in the construction of b; computations for c are quite 
similar. _ 

Example 1. To find b we can bind (w(J(t) + 5 J, #) - |y ( J(r)) u(j(r) + 6 J, tf)) 2 , 

which has the form /i(J(r), <5J, i?) with ft, a polynomial as above; the square root of 
the majorant fc(j(r),r) is b{r,r). 

Example 2. This computation is very similar to the one for Example 1. 
Example 3. The left hand side of Eq. (2.24) is 

H J(r) + 8 J, 0) - 2( J(r))v( J(r) + 5 J, 0) | = (E.2) 



3-4cos^ + cos(2^) < 2 < 



=|W| 



4(j(r) + 5J) 3 ^ (J(r) + 5J) 3 (J(r) - r) 3 
Example 4. In this case, 

|«( J(r) + tJ, *) - f ( J(r)) ,( J(r) + W, *)|' = (J , (T)+ ^ ( "ff +fja)< (E.3) 

with h a polynomial as before. After finding for h a bound of the form (E.l), we 
combine it with the obvious relation (J*(t) + 5J*) -4 ^ (J 1 (t) — ^)~ 4 for r := |<5J|; 
the square root of the final majorant is 6(r, r). 

iii) Up to now, we have not explained how to get elementary bounds of the form 
(E.l) on a polynomial h. Here we illustrate a general procedure (computations to 
apply it in Examples 1-4 are generally too tedious to be made by hand, but are 
easily implemented on MATHEMATICA). 

iii a) In the expression of h(J, 5 J, •&), if d — 1 we put 5 J = r cosip with ijj = or n; 
if d — 2, we set 5 J = (r cos ip,r sin ip) with ip e T. 

iii b) Now, ft(J, <5J, i?) has the form of a trigonometric polynomial in ^ with 
coefficients depending on J. We write this in a canonical form, reexpressing any 
term in d and ip as a linear combination of sines and cosines (e.g., cos 4 1? sin(2i9) 2 = 
(l/32)(5 + 4cos(2tf) -4cos(4^) -4cos(6^) - cos(8tf)). 

iii c) As a final step, we bind each summand of h using the relations | cos |, | sin | ^ 1. 
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